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ABSTRACT 

We present 74 MHz radio continuum observations of the Galactic center re- 
gion. These measurements show nonthermal radio emission arising from molec- 
ular clouds that is unaffected by free-free absorption along the line of sight. We 
focus on one cloud, GO. 13-0. 13, representative of the population of molecular 
clouds that are spatially correlated with steep spectrum (a|27MHz = 1-3 ± 0.3) 
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nonthermal emission from the Galactic center region. This cloud lies adjacent to 
the nonthermal radio filaments of the Arc near 1~ 0.2° and is a strong source of 
74 MHz continuum, SiO (2-1) and Fel Ka 6.4 keV line emission. This three-way 
correlation provides the most compelling evidence yet that relativistic electrons, 
here traced by 74 MHz emission, are physically associated with the GO. 13-0. 13 
molecular cloud and that low energy cosmic ray electrons are responsible for the 
Fel Ka line emission. The high cosmic ray ionization rate ~ 10~ 13 s _1 H^ 1 is 
responsible for heating the molecular gas to high temperatures and allows the 
disturbed gas to maintain a high velocity dispersion. LVG modeling of multi- 
transition SiO observations of this cloud implies H 2 densities ~ 10 4-5 cm -3 and 
high temperatures. The lower limit to the temperature of GO. 13-0. 13 is ~ 100K, 
whereas the upper limit is as high as 1000K. Lastly, we used a time-dependent 
chemical model in which cosmic rays drive the chemistry of the gas to investigate 
for molecular line diagnostics of cosmic ray heating. When the cloud reaches 
chemical equilibrium, the abundance ratios of HCN/HNC and N 2 H + /HCO + are 
consistent with measured values. In addition, significant abundance of SiO is 
predicted in the cosmic ray dominated region of the Galactic center. We dis- 
cuss different possibilities to account for the origin of widespread SiO emission 
detected from Galactic center molecular clouds. 

Subject headings: astrochemistry - ISM: abundances — ISM: molecules — ISM: 
lines and bands 

1. Introduction 

The interstellar medium (ISM) of the nucleus of our Galaxy has unique characteristics. 
Gas clouds are subject to the strong tidal field of the nucleus and are required to have higher 
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density than in the Galactic disk in order to become gravitationally unstable and form stars. 
This region of the Galaxy hosts two components in its central few hundred parsec. One 
is the reservoir of molecular gas showing enhanced molecular emission with higher velocity 
dispersion (~ 20 — 30) km s" 1 and gas temperature (50-200K) compared to elsewhere in the 
rest of the Galaxy 1 " 7 . The other is the prevalence of a mixture of thermal emission from 
ionized gas at a temperature of few thousand degrees and nonthermal radio synchrotrone 
emission from relativistic electrons (e.g.; Nord et al. 2004; Law et al. 2008). The synchrotron 
emission is best viewed at low radio frequencies (< 1GHz) whereas thermal free-free emission 
is better detected at high frequencies (> 1GHz). Large-scale surveys of this region suggest 
that relativistic electrons and molecular gas co-exist and possibly interact with each other. 
However, the apparent correlation could result from chance coincidence along the long line of 
sight towards the Galactic center. We present the first evidence for low frequency nonthermal 
emission closely tracing molecular gas, establishing that cosmic ray electrons are physically 
associated with individual molecular clouds. We study one molecular cloud GO. 13-0. 13 in 
detail and postpone the discussion of other Galactic center clouds to elsewhere. 

The interaction between cosmic-ray particles and molecular gas has several far reaching 
consequences. For one, cosmic rays play an important role in star formation processes as 
they are the primary source of ionization in molecular clouds that are self-shielded from UV 
radiation field. The interaction of cosmic ray electrons heats the gas to a higher temperature, 
which increases the Jeans mass and causes the initial mass function (IMF) to become top- 
heavy 8 " 10 . The higher ionization fraction due to the impact of these electrons reduces the 
damping of MHD waves and helps to maintains the high velocity dispersion of molecular 
gas in the nuclear disk. In addition, this interaction strengthens the coupling of gas to the 
magnetic field slowing star formation by increasing the time scale for ambipolar diffusion 
before the onset of gravitational collapse. 
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Recent measurements indicate a vast amount of H3 and H 3 + 11-14 as well as high 
temperature molecular gas distributed in the Galactic center. The inferred minimum cosmic 
ray ionization rate, ( ~ 10~ 15 s _1 H -1 , is one to two orders of magnitude higher in the 
Galactic center region than in the Galactic disk 10,13 . The elevated cosmic ray ionization 
rate will increase the ionization fraction of electrons in molecular gas, and drive ion-neutral 
chemistry. Another consequence is enhanced fluorescent Fel Ka 6.4 keV emission 15 ' 16 . This 
nonthermal X-ray emission is uniquely detected from Galactic center clouds and results from 
filling the inner K-shell vacancies of neutral iron created by the impact of low energy cosmic 
ray electrons or by irradiation by X-ray photons 17-20 . A widely accepted model argues that 
the 6.4 keV line is a result of irradiation by a hypothetical transient source associated with 
Sgr A*, which was active about 400 and again 100 years ago, and that we are seeing this 
transient's light echo in the 6.4 keV emission 21 . We discuss another scenario in which low 
energy cosmic rays can contribute significantly in production of the Fel Ka line emission 
from Galactic center clouds. 

The Galactic center population of molecular clouds have physical conditions similar to 
hot cores, but with the exception of Sgr B2, there is no strong evidence for ongoing massive 
star formation throughout the inner few hundred parsecs of the Galaxy, a region known as 
the central molecular zone (or CMZ) 22 . Line emission from a variety of molecular species 
has been detected throughout the CMZ including those that are produced in the gas phase 
(e.g.; HCN, HCO + , HNC, N 2 H + ) as well as those processed on grain surfaces, e.g. SiO 6,r . 
With the exception of the strong radiation field near Sgr A* and the Arches cluster, PDRs 
can not be very important in the dense self-shielded clouds. In addition, the X-ray flux from 
the Galactic center is too weak to qualify for the application of XDR models. Global heating 
by cosmic rays can explain why the gas temperature is significantly higher than the dust 
temperature in a large fraction of the gas and dust clouds in the nuclear disk. In the cosmic- 
ray dominated region of the Galactic center, cosmic-ray heating should be more significant 
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than ambipolar and turbulent heating 23 . Thus, it is natural to examine the chemistry of the 
gas in the cosmic-ray dominated region and attempt to identify molecular line diagnostics 
of cosmic ray heating. 

To this end, we present high resolution 74 MHz observations of the Galactic center 
and show the spatial correlation between molecular gas and 100 MeV cosmic ray electrons. 
Enhanced 74 MHz emission appears to coincide with a large concentration of molecular 
gas in the inner 5° of the Galactic center 24 . Here, we focus on GO. 13-0. 13, a representative 
molecular cloud, and compare its 74 MHz emission with that of excited rotational transitions 
of SiO and CS molecules. This molecular cloud lies along the nonthermal filaments of the 
radio Arc near 1~ 0.2°, the most prominent network of magnetized filaments in the Galactic 
center. Molecular line emission from CO, CS, SiO, H 13 CO + and CS suggests that this cloud 
has a high column density and gas temperature 25-27 . The kinematics of CS line emission 
from GO. 13-0. 13 suggests an expansion of molecular gas into the nonthermal filaments 25 . CO 
observations imply gas temperature T> 70K and column density N(H 2 )=6-7x 10 23 cm -2 . 
Multiple transitions of NH3 line emission from this cloud have measured two temperature 
components giving a range T~25K and T~ 125 — 200K 28 . The low temperature component 
of molecular gas is similar to the measured 18-22 K dust temperature in clouds in the inner 
2° x 1° of the Galaxy 29 ' 30 . 

The structure of this paper is as follows. We first show 74 MHz emission from Galactic 
center molecular cloud GO. 13-0. 13 and estimate the variation of spectral index of nonthermal 
emission between 327 and 74 MHz. We then show the association of the nonthermal radio 
emission from GO. 13-0.13 with Fel Ka line emission at 6.4 keV and estimate the cosmic 
ray ionization rate needed to produce the Fel Ka line emission. We will also present new 
molecular line observations of four rotational transitions of SiO to determine the temperature 
and density of molecular gas in GO. 13-0. 13 more accurately. We consider the interaction 
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of relativistic electrons with the gas and compute the cosmic ray ionization rate, implied 
by the 74 MHz emission, as a function of assumed magnetic field strength. We present 
the dependence of the total cooling rate of the gas on gas temperature for high molecular 
gas densities. Lastly, we study the chemical consequences of the interaction of cosmic ray 
electrons with the molecular gas by modeling the abundance ratios of several molecular 
species as a function of time. A time- dependent gas-grain chemical model 31 is used to 
explore how abundance ratios of five representative molecular species vary with gas density 
and cosmic-ray ionization rate. 

Our study of GO. 13-0. 13 is not only relevant to understanding the molecular component 
of the nuclear disk of our Galaxy 8 but also to external galaxies and ultraluminous infrared 
galaxies (ULIRGs) where cosmic rays are thought to be the driving mechanism for star 
formation 9 . 

2. Observations and Data Reduction 

We imaged the Galactic center at 74 MHz and mapped molecular line emission from 
several molecular clouds in the Galactic center using Mopra, CSO and CARMA. We only 
present observations of GO. 13-0. 13 here and will describe details of observations of other 
Galactic center clouds elsewhere. 

2.1. VLA 

The 74 MHz radio continuum observations of the Galactic center were described previously 24 . 
Briefly, these measurements were taken in multiple configurations of the Very Large Ar- 
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ray (VLA) of the National Radio Astronomy Observatory^] (NRAO), with a resolution of 
114" x 60" (PA= —5°) and rms sensitivity of ~ 0.12 Jy. These observations are unprece- 
dented in spatial resolution and sensitivity at these low frequencies. Given the strong fre- 
quency dependence (z^ 2 ' 1 ) of free-free absorption, foreground and embedded thermal sources 
are seen in absorption against the strong nonthermal emission from the Galactic center. The 
suppression of thermal emission due to high opacity of ionized gas at 74 MHz allows us to 
readily identify nonthermal sources that emit at this frequency. 



2.2. CSO 

Observations of the SiO (5-4) and (6-5) rotational transitions at 217.1 and 260.5 GHz, 
respectively, were carried out in 2010 September using the wideband 230 GHz "Z-Rx" re- 
ceiver of the Caltech Submillimeter Observatory (CSO) on Mauna Kea, Hawaii. The weather 
conditions were good, characterized by a 225 GHz zenith opacity of ~ 0.05 — 0.07, equivalent 
to 1—1.5 mm of precipitable water vapor. Typical single sideband system temperatures at 
the relatively low elevation of the Galactic center were ~250 K. The CSO FWHM beam 
size at the two frequencies is ~ 34" and 28", respectively, and the main beam efficiency, 
as determined from total power observations of Jupiter, was ~ 70%. Spectra were taken 
in the "on-the-fly mapping" mode, on a ~ 15" grid, using a designated off position at 
(a, 5)2ooo.o =(17 /l 45 m 13. s 2, -28°45 / 31."0), out of the Galactic plane. Multiple maps scanned 
in orthogonal directions were averaged together to avoid striping. The final per pixel in- 
tegration time (ON) was 14-17 sec for SiO 5-4 and 25-31 sec for SiO 6-5. We used the 
high-resolution CSO facility FFTS spectrometer with 8192 channels and a total bandwidth 
of 1 GHz. 

1 The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated 
under a cooperative agreement by Associated Universities, Inc. 
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2.3. Mopra 

We observed GO. 13-0. 13 on 2010, June 28 using Mopraj^] in the on-the-fly mapping 
mode. The map size is 5' by 5', and is centered at Galactic longitude and latitude (1, 
b) = (0.1092°, —0.1000°). We used a latitude scan direction (orthogonal to the Galactic plane), 
and a row spacing of ~ 15". We used the MOPS spectrometer in Wideband mode, with a 
central frequency of 89.690 GHz. This results in a frequency coverage of ~ 85.6 — 93.8GHz 
across four sub-bands. Each sub-band has 8192 channels, a channel width of 269.5 kHz 
(~ 0.9 km s _1 ), and two linear polarizations. The final map is the average of the two 
polarizations. Regular pointing checks done throughout the observations, and were within 
~ 4" (the beam size is ~ 38" at these frequencies). The intensity is given in T^ which can 
be divided by the main beam efficiency 0.49 to convert the intensity to T m &. 

We also made deep pointed observations toward four positions (A-D) within GO. 13- 
0.13 on 2010, June 30 (see Table 1). These spectra were obtained in an on-off position 
switching mode for a total on-source time of 288 s. The MOPS spectrometer was configured 
in wideband mode, identically to the setup of the map described above. A total of 16 
molecular lines CH 3 CN, C 13 CN, 13 CS, N 2 H+, HNC, HCCN, HNCO, HCN, HCO+, 28 SiO, 
HC 15 N, SO, H 13 CN, H 13 CO, SiO, NH 15 C and C 2 H are detected from four positions in GO. 13- 
0.13 (see Fig. 7). 

2 The data was obtained using the Mopra radio telescope, a part of the Australia Telescope National 
Facility which is funded by the Commonwealth of Australia for operation as a National Facility managed by 
CSIRO. The University of New South Wales (UNSW) digital filter bank (the UNSW-MOPS) used for the 
observations with Mopra was provided with support from the Australian Research Council (ARC), UNSW, 
Sydney and Monash Universities, as well as the CSIRO. 
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2.4. BIMA 

We mapped GO. 13-0. 13 in CS(2-1) at 97.980968 GHz and HCO+(1-0) at 89.188518 
GHz with the Berkeley-Illinois-Maryland Association (BIMA) interferometer]^] during the 
2002-2003 observing season. The BIMA array was a mm-wavelength interferometer located 
in Hat Creek, California 32 , consisting of ten 6.1 m antennas. These antennas have since 
been combined with the Owens Valley Radio Observatory 10 m antennas into the CARMAP] 
millimeter array at Cedar Flat, California. 

The cloud was observed with 20 mosaicked pointings on a hexagonal grid corresponding 
to Nyquist sampling of the primary beam and covering a total area of about 5' x 6'. The 
pointing center was (a, 5)2000.0 =(17 h 46 m 24.'0, -28°45'00."0) with a velocity of 30 km s" 1 
with respect to the local standard of rest. The HCO + and CS data were obtained simul- 
taneously in opposite sidebands. Each of the spectral lines were observed with a 50.0 MHz 
correlator window with 256 channels. Along with the spectral data, continuum windows 
with a total coverage of about 800 MHz were also observed. We observed in the C and 
B configuration of the array, sampling spatial frequencies from 1 to 78 kA. Absolute flux 
calibration was derived from observations of Mars or Uranus immediately before or after 
the source track. The quasar 1733-130 was used as the phase calibrator and secondary flux 
calibrator. To fill the inner uv spacings, we also obtained fully-sampled maps covering the 
region imaged with the array in the HCO + , and CS using the Five Colleges Radio Astronomy 
Observatory 14 m telescope, as well as the previously published CS(2-1) maps made with 
Nobeyama 45-m radio telescope (NRO) 25 . A more detailed high resolution kinematic study 



3 The BIMA interferometer was operated under a joint agreement between the University of California, 
Berkeley, the University of Illinois, and the University of Maryland with support from the National Science 
Foundation. 
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Combined Array for Research in Millimeter-wave Astronomy; 



http://www.mmarray.org 
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of GO. 13-0. 13 including the results of HCO + observations will be given elsewhere. 

3. Results 

Figure la,b shows the distribution of 74 MHz continuum and the CO (3-2) line from 
the inner 1.5° x 12' of the Galactic center. A layer of molecular gas in the CMZ runs parallel 
to the Galactic plane and is seen in the distribution of CO (3-2) line emission 33 . This layer 
includes some of the most prominent molecular clouds in the CMZ such as Sgr B2 and Sgr 
Bl, the 40, 50, -30 and 20 km s _1 Sgr A and Sgr C clouds. Some of the prominent clouds 
are labeled on Figures la,b. At positive longitudes, the Sgr B2 complex, GO. 36-0. 09 and 
GO. 13-0. 13 have 74 MHz counterparts and at negative longitudes the layer of molecular 
gas that runs parallel to the Galactic plane G359.75-0.23 shows a 74 MHz counterpart. 
Another presentation of the CO (3-2) and 74 MHz emission is shown in Figure lc. The 
distribution of 74 MHz emission appears to follow the CO (3-2) line emission except in two 
types of regions. One type is associated with holes created by free-free absorption suppressing 
background 74 MHz emission and the other is where nonthermal continuum sources have 
no molecular counterparts. Overall, the 74 MHz emission shows similar morphology to 
that of CO (3-2) line emission apart from absorption features manifesting as holes in the 
74 MHz distribution. To clarify the correlation between absorption features and free-free 
emission, Figure Id presents a continuum map at 8.5 GHz based on GBT observations 34 . 
Prominent continuum sources are labeled on this figure. Holes in the 74 MHz continuum map 
coincide with strong continuum emission at 8.5 GHz, consistent with the free-free absorption 
coefficient increasing significantly at low frequencies. Thus, thermal sources become optically 
thick and are seen in absorption against strong nonthermal background emission 35 whereas 
nonthermal sources that are not contaminated by significant foreground thermal emission 
are clearly identified at 74 MHz. The morphological comparison of these figures suggest that 
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nonthermal continuum emission at 74 MHz arises from Galactic center molecular clouds, 
implying that cosmic-ray electrons are interacting with molecular gas in this unique region 
of the Galaxy. 

We now focus on the GO. 13-0. 13 cloud. Figure 2a shows contours of CS (1-0) molecular 
line emission from this cloud superimposed on a grayscale 5 GHz continuum image. A 
network of linearly polarized filaments of the radio Arc running almost perpendicular to the 
Galactic plane at a PA~ —170° lies at the eastern edge of the cloud 37 . The morphology of 
molecular line and radio continuum emission from GO. 13-0. 13 suggests that the molecular gas 
is surrounded by, and is interacting with, nonthermal radio filaments. The kinematics of the 
molecular gas suggests that expansion of the cloud is responsible for a dynamical interaction 
between the eastern limb of the GO. 13-0. 13 and the magnetized nonthermal filaments 25 . 
Figure 2b shows high resolution CS (2-1) line emission from the GO. 13-0. 13 molecular cloud 
(red) superimposed on a 1.4 GHz continuum image (green). The molecular gas distribution is 
clumpy with the morphology of a "boot" , and is edge-brightened parallel to the nonthermal 
radio filaments of the Arc near 1~ 0.2° on the eastern edge of the cloud. We note that the 
western edge of the cloud lies parallel to another magnetized nonthermal filament GO. 087- 
0.087 38 . The molecular line images show a ridge of emission that deviated from a straight 
line perpendicular to the Galactic plane (labelled as "Meanderning feature" on Figure 2a) 
near 1~ 8', b~ —7'. This ridge of emission appears to curve around a nonthermal filament 
near 1= 8' 4", b=— 6' 55". This high resolution map of molecular line emission provides 
the strongest morphological support for the interaction of the nonthermal filaments and the 
eastern edge of GO. 13-0. 13. 

Figure 3a,b show contours of 74 and 327 MHz emission, respectively, superimposed on a 
grayscale image at 5 GHz. The 74 MHz emission is produced mainly by nonthermal processes 
whereas 5 GHz and 327 MHz emission can trace emission produced by both thermal and 
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nonthermal processes. The distribution of 74 MHz emission differs remarkably from the 327 
MHz and 5 GHz continuum. The 74 MHz emission arises from the vertical filaments as well 
as the region where the GO. 13-0. 13 molecular cloud lies. There is no evidence of thermal 
emission from GO. 13-0. 13 at 327 MHz or at higher frequencies, thus, the spectrum derived 
here is not contaminated by thermal emission. The vertical filaments seen in the grayscale 
image at 5 GHz have counterparts at 74 MHz except in the thermal region to the north of 
the filaments near GO. 18-0.04. This is because GO. 18-0.04 is an HII region dominated by 
free-free emission at high frequencies, thus the lack of 74 MHz emission is due to significant 
free-free absorption 35 . 

What is remarkable about images shown in Figure 3a,b is the discovery that enhanced 
nonthermal emission at 74 MHz arises from the molecular cloud GO. 13-0. 13. The spectrum 
of emission from the center of the cloud must be steep given that there is no strong 327 
MHz emission from GO. 13-0. 13. Figure 3c shows the spectral index distribution (where the 
flux density F„ oc i*~ a ) between 327 and 74 MHz. The spectrum is inverted (a < 0) at 
the position of the filaments and is steep (a ~ 1.3) where the molecular line and 74 MHz 
emission peak. There is no evidence of thermal emission from GO. 13-0. 13 at 327 MHz or at 
higher frequencies, thus, the spectrum derived here is not contaminated by thermal emission. 

Another version of the spectral index distribution is made by comparing the background 
subtracted profiles of the emission at 74 and 327 MHz. Figures 3d,e show two different cross 
cuts made at constant longitude 1=— 6' and constant latitude b=— 8', respectively. This 
technique follows earlier work by Law et al. (2008). The spectral index distribution, as 
shown in Figure 3d, is consistent with a steep spectrum (a = 1.1 ± 0.2) becoming flatter or 
inverted at the nonthermal vertical filaments. Figure 3e shows a « —0.2 ± 0.05 between 74 
and 327 MHz at the position of nonthermal filaments. The electron energy spectrum of the 
center of GO. 13-0. 13 is remarkably steep, p^3.2 (where p=2a+l corresponding to energy 
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spectrum E p ). In contrast, the spectral index of the vertical filaments is remarkably flat 
pf»0.6. The variation of the spectral index of nonthermal emission from the filaments and 
GO.13-0.13, is Aa ~ 1.5. 

GO. 13-0. 13 is one of the Galactic center molecular clouds showing bright Fel Ka line 
emission at 6.4 keV 39 . Figure 4 shows contours of fluorescent Fel Ka line emission based 
on Chandra observations convolved to a resolution of 30" and superimposed on a grayscale 
distribution of integrated CS (1-0) line emission. High resolution Chandra and CARMA 
images of this cloud show a clumpy distribution of Ka line emission 20 tracing CS (1-0) line 
emission. The overall distribution of Fel Ka line emission arising from GO.13-0.13 is similar 
to that of low frequency radio emission molecular line CS and SiO emission. The distribution 
of the Ka line emission is clumpy and shows a ridge of emission adjacent to a nonthermal 
feature at 1=8' 10" b=— 6' 55" . Remarkably, this X-ray ridge feature is similar to that of 
CS emission (cf. Meandering feature in Fig. 2b). Both the X-ray and CS ridges curve 
around a nonthermal filament. The morphology of 6.4 keV line emission supports the idea 
that nonthermal radio filaments are interacting with the edge brightened Fel Ka emitting 
molecular cloud GO.13-0.13 at 6.4 keV. These images suggest a three-way correlation between 
cosmic ray electrons, molecular gas and Fel Ka line emission. 

In order to estimate the gas density and temperature in GO.13-0.13, we measured the 
line intensities of four transitions of SiO (6-5), (5-4), (2-1) and (1-0). We used the SiO 
(1-0) and (2-1) images taken from Handa et al. (2006) and convolved them to the same 
resolution and velocity coverage, between -10 to 50 km s^ 1 as those of SiO (6-5) and (5-4) 
transitions. Figure 5a shows contours of velocity integrated SiO (5-4) line emission mapped 
over a limited area of GO.13-0.13 superimposed on a 5 GHz continuum image. We applied 
a Large Velocity Gradient (LVG) model to the two E and NW rectangular boxes, as shown 
schematically, in order to derive the physical parameters of molecular gas in GO.13-0.13. 
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The LVG grid was computed for a constant SiO column density of 1.6xl0 13 cm~ 2 and a 
line width of 30 km s _1 . This approximately reproduces the observed line intensities for the 
range of densities and temperatures applicable to the gas in the Galactic center region. With 
these assumptions, the emission is optically thin (typical optical depths ranging between 0.1 
and 10~ 3 , depending on the transition). Thus, optical depth effects do not influence the 
LVG model results as long as the beam filling factor is close to unity. A smaller beam filling 
factor will change the optically thin assumption we have made. Figure 5b, c show the density 
and kinetic temperature constraints estimated from three SiO line ratios of the E and NW 
boxes in three different colors, respectively. The three line ratios give a wider range of gas 
temperature between ps 100 — 1000K for box E and and « 400 — 1000K for box NW. A lower 
range of gas temperatures > 100K is found from the region near the nonthermal filaments 
of the Arc, box E as shown in Figure 5a, than that of the peak emission in the NW box 
having T>400K. The implied kinetic temperatures are consistent with earlier measurements 
based on multiple transitions of NH 3 28 . Additional constraints have also been made by 
estimating molecular gas pressure from the NW box of Figure 5a. The top three panels of 
Figure 5d show the plots of line ratios for all models in the grid as function of logarithm of 
pressure (nT). The bottom panel shows a sum of [(model-observed) /observed] 2 for the three 
line ratios and the the best fit is obtained for log(nT) in the range of 6.4-6.9. The line ratios 
give values of density and temperature that are consistent with log(n)~ 4.4 and log(T)~ 2.3. 
for minimum \ 2 values. Similar values were obtained for the E box of Figure 5a. 

Tsuboi et al. (2011) used the SiO (2-l)/(l-0) line ratios to estimate a molecular gas 
density of ~ 10 4 cm -3 at a kinetic temperature of 60K. On the other hand, Oka et al. (2001) 
used the CO (2-1) / (1-0) intensity ratios and argued that GO. 13-0. 13 has a low density 
~ 10 2 cm~ 3 , high kinetic temperature > 70K component. Overall, our model fits which 
uses multiple transitions of SiO are suited to measure gas density and temperature, though 
with systematic uncertainties from using different telescopes. The application of the LVG 
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excitation code constrains the gas density of hydrogen nuclei n H ~ (1 — 3) x 10 4 cm -3 and 
temperature T~ (1 — 2) x 100K but we can not rule out gas temperature as much as 1000K 
and lower densities n~ 10 3 cm" 3 . 

Two studies have made discrepant estimates of the column density of the gas in GO. 13- 
0.13 26 ' 40 Using kinetic temperature T fc =70K, Handa et al. (2006) used H 13 CO+ (J=l-0) 
line emission which is expected to be optically thin with abundance of 10" 10 to derive an 
H 2 column density (6-7) xlO 23 cm" 3 . On the other hand, Amo-Baladron et al. (2009) 
determined the physical condition by applying an LVG model to the SiO (2-1) and (3-2) 
transitions and then estimated the column of H 13 CO + (J=l-0). Using the same abundance 
as that of Handa et al. (2006), the derived column density of molecular hydrogen from this 
study is (2-8) xlO 22 cm" 3 . To resolve this discrepancy, we estimated the column density of 
molecular gas by using submillimeter emission from dust grains in GO. 13-0. 13. We used the 
data taken from Pierce-Price et al. (2000) who derived 513 M Jy -1 beam -1 at 850yum 
with the assumption that dust temperature is 20K, metallicity twice solar and opacity index 
(3=2. Using background subtracted flux densities of 0.5 and 1.5 Jy from the NW and E 
boxes at 850/iin and a beam of 15", we find the column density of molecular hydrogen ranges 
between ~ 3.5 x 10 22 and ~ 10 23 , cm -2 respectively. There is uncertainty in accurately 
determining the background emission. Given this uncertainty, our estimate of the column 
density falls between previous measurements of Handa et al. (2006) and Amo-Baladron et 
al. (2009). Herschel maps presented by Molinari et al. (2011) give column density values 
that are consistent with our estimates. 



4. Discussion 



Our interest in studying Galactic center molecular clouds stems from three earlier stud- 
ies, all of which had suggested the interaction of cosmic ray electrons with molecular gas and 
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that cosmic ray electrons are enhanced in this region. First, fluorescent 6.4 keV Ka iron line 
emission from several molecular clouds in the Galactic center may be the consequence of the 
interaction of low energy cosmic ray electrons with molecular gas 20 ' 41 , or alternatively, the 
irradiation of molecular clouds by a hypothetical X-ray flash associated with past activity 
of Sgr A* 21 . Second, nonthermal bremsstrahlung from this electron population can also ex- 
plain the diffuse 7-ray emission from the central 2° x 1° of the Galactic center 42 . Third, the 
strong H3 absorption along several lines of sight towards the Galactic center imply that the 
ionization rate is at least an order of magnitude higher than elsewhere in the Galaxy 10 . This 
is consistent with the low energy cosmic ray electron production of the Fel Ka line emission. 
Finally, low frequency diffuse radio emission 43 combined with spectral index measurements 23 
suggest strong nonthermal continuum emission from the Galactic center. Our detailed study 
provides the most compelling evidence yet that relativistic electrons, here traced by 74 MHz 
emission, are physically associated with the GO. 13-0. 13 molecular cloud. This cloud is one 
of the strongest sources of Fel Ka line emission in the Galactic center. We complete our 
analysis by studying the complex morphology of GO. 13-0. 13, estimating the cosmic ray ion- 
ization rate, the expected Fel Ka line emission at 6.4 keV, and modeling chemical signature 
of the interaction. 

4.1. Morphology & Kinematics 

GO. 13-0. 13 is a complex molecular cloud with kinematics that suggest expansion into 
the vertical filaments of the radio Arc 25 Morphologically, this cloud lies at the center of a 
circular-shaped structure, known as the radio Arc bubble, apparent in mid-IR images taken 
with the Spitzer and MSX observatories 44 " 47 . Figure 6 shows contours of SiO (2-1) line 
emission, obtained with Mopra, superimposed on a three-color Spitzer IRAC image. The 
Arc bubble is noted exterior to the south of GO. 13-0. 13. The origin of the bubble and its 
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relation to the radio Arc and the GO. 13-0. 13 cloud is not well understood. Past studies have 
suggested that the radio Arc bubble lies in the vicinity of the Quintuplet star cluster and is 
produced by stellar winds or supernova explosions sweeping up interstellar material 47 . The 
nonthermal radio emission at 74 MHz and expansion of GO. 13-0. 13 point to the possibility 
that the Arc bubble may be produced by the same event. It is possible that an energetic event 
is driving a shock into GO. 13-0. 13. The shock reaches the edge of the cloud and before it 
encounters much lower density of the gas exterior to GO. 13-0. 13, then sweeps through the low 
density gas and creates a edge-brightened bubble. In this picture, the expansion of GO. 13- 
0.13 into the nonthermal vertical filaments accelerates particles along the filaments 25,48 , thus 
generating a young population of electrons running along the magnetized filaments with an 
unusually flat energy spectrum (p~0.6). The puzzle, however, is the origin of the steep 
energy spectrum of particles with p> 3 in GO. 13-0. 13. This is because energy losses are 
increasingly severe for lower energy electrons, tending to flatten their energy spectrum at 
low energies. 



4.2. Cosmic Ray Ionization Rate 

Using equation (3) of Yusef-Zadeh et al. (2012), the cosmic ray ionization per second 
per hydrogen nucleus, (, depends on the observed synchrotron intensity l u , the magnetic field 
B, the depth of the source of emission along the line of sight, L, and the energy spectrum of 
the electrons with the index p, assumed to be power-law E~ p between 0.1 MeV and 1 GeV 
and a — {p — l)/2 is the synchrotron spectral index and l v oc v~ a . 

„ 3.1 x IP-* 4 l v fJL1 - - 1 / B B -i h-1 m 

S ~ p-1 Jyarcmin" 2 V GHz/ VSOpc; V 100 /^/ ' 



Figure 7a shows the inferred cosmic ray ionization rate as a function of the magnetic 
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field strength for different values of a. We assumed that L~2.5 pc (1' corresponds to 2.4pc 
at the Galactic center distance 8.5 kpc) and the observed surface brightness of 0.71 Jy 
arcmin -2 at 74 MHz averaged over the inner 5' x 5' of GO. 13-0. 13. The black dot on the 
curve gives the equipartition magnetic field for the electron population producing diffuse 
synchrotron emission in the dense molecular cloud GO. 13-0. 13. The ionization rate due to 
electrons increases for steeper spectral index values for a fixed emissivity at the observed 
frequency, as there is a successively larger population of lower-energy electrons radiating at 
lower frequencies. The equipartition magnetic field ranges between 30 /zG and 0.3 mG for a 
values between 0.25 to 1.5. The corresponding ( is estimated to be ~ 4 x 10~ 15 and 10~ n 
s -1 , respectively. A cosmic ray ionization rate > 10~ 13 s _1 H _1 is sufficient to fully dissociate 
the gas on a time scale of < 10 6 years. To avoid this, the magnetic field has to be larger 
than lmG for a = 1.25. 

4.3. The 6.4 keV Neutral Iron Line Emission 

We apply the cosmic ray model to the GO. 13-0. 13 molecular cloud by estimating the 
Ka line emission from the interaction of the low energy nonthermal electrons responsible for 
synchrotron radio emission detected at 74 MHz. Using equation (6) of Yusef-Zadeh et al. 
(2012), 

I Ka « 8.7 x 10- 8 (y^tt) (1^=2) P h ^ cm " 2 arcmin " 2 > ( 2 ) 

for N(H2) = (1 — 3) x 10 23 cm -2 , then the cosmic ray ionization rate needed to give the 
background subtracted Ka line flux of 3 x 10~ 6 ph cm~ 2 s _1 arcmin -2 derived from Suzaku 
observations 39 is ( ~ (1 — 3) x 10~ 13 s -1 , respectively. A value of cosmic ray ionization rate 
( ~ 10~ 13 s" 1 is indicated by a dashed line on Figure 7a. Using the spectral index values 
between a = 1 and 1.25 and the gas density between few times 10 4-6 cm -3 , the magnetic 
field strength is constrained to values between ~ 0.5 and 1 mG, respectively. 
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A widely accepted alternative model argues that the 6.4 keV Fel Ka line emission results 
from irradiation of molecular clouds near the Galactic center by a hypothetical transient 
source associated with Sgr A*, which was active about 400 and again 100 years ago 21 . The 
variability of 6.4 keV line emission on a short time scale provided a strong evidence in 
support of this hypothesis. However, the evidence for nonthermal emission from molecular 
clouds, as presented here, also predicts time variability in the context of the cosmic ray 
picture 42 . A supernova event expanding inside the cloud could change not only the ionization 
characteristics of the molecular gas in GO. 13-0. 13 but also causes the observed 6.4 keV time 
variability. Thus, the time variability of Ka line emission need not be clear support for the 
X-ray irradiation model. 



4.4. Cosmic Ray Heating of Molecular Gas 

To estimate the temperature of molecular gas subject to high levels of cosmic-ray ioniza- 
tion, we extended our previous estimates of the cooling rate 42 to higher densities, using the 
fitting formula for cooling by rotational transitions of CO, H 2 , and H 2 provided by Neufeld 
& Kaufman (1993) and Neufeld, Lepp & Melnick (1995). We also included an LVG calcula- 
tion of the contribution by OI fine-structure lines, using the atomic data and collision rates 
in the Leiden Atomic and Molecular Database 49 . We adopted an abundance of CO relative 
to H 2 of 2.8 x 10~ 4 , twice the values adopted by Neufeld, Lepp & Melnick (1995) to reflect 
the higher metallicity of the gas in the CMZ. The abundance of OI and H 2 is temperature- 
dependent because of the conversion of OI and 2 to H 2 by neutral-neutral reactions at 
temperatures above 200 K, so we assume twice the chemical equilibrium abundances of OI 
and H 2 found by Neufeld, Lepp & Melnick (1995): OI/H 2 = 2 x 1(T 4 /(1 + (T/220K) 14 ), 
and H 2 0/H 2 = dex(-6.222 + 2.831/(1 + (245K/T) 14 )). We adopted a cloud column of 
2 x 10 23 H 2 cm- 2 and a line FWHM 20kms- 1 . 



- 20 - 



The resulting total cooling rate per H 2 molecule is plotted in Figure 7b for representative 
H 2 densities of 10 3 , 10 4 and 10 5 cm~ 3 . The cooling rates on the left hand axis are mapped 
to the right-hand axis showing the cosmic ray ionization rate that would supply heat at 
the same rate, assuming that each ionization is associated with the deposition of 12.4 eV of 
heat 42 ' 52 . The dip in the cooling rate at T 250 K is due to the increasing importance of 
01 and H2O cooling at higher densities and the switch from 01 to H2O at this temperature. 
We conclude that an ionization rate of 10~ 14 to 10 _13 s _1 would yield gas temperatures in 
the range 50-200 K, the upper range consistent with the temperatures estimated using LVG 
shown in Figure 5b, c. 



4.5. High Velocity Dispersion of Molecular Clouds 

High cosmic-ray fluxes in molecular clouds affect star formation by heating the gas and 
increasing its ionization fraction. Higher cloud temperatures increase the Jeans mass, po- 
tentially changing the IMF, while high ionization increases magnetic coupling to the cloud 
material, reducing ambipolar diffusion and increasing the time taken for gravitationally un- 
stable cores to contract to the point that they overwhelm their magnetic support. Another 
consequence of increased ionization is the therefore reduced damping of MHD waves, con- 
tributing to sustaining Alfvenic velocity fields within the clouds, which may assist in explain- 
ing the observed high velocity dispersion of molecular clouds in the nuclear disk (e.g. Oka et 
al. 1998; Martin et al. 2004). In a weakly ionized medium, waves with frequencies u ~ kvA 
below the collision frequency of neutral particles with ions, u ni = rii < av >, are damped on 
a time scale 2v ni /uj 2 53 ' 54 directly proportional to rii which is increased by a factor of ~ 100 
compared to the Galactic disk for ( ~ 10~ 13 s -1 . The power input required to maintain wave 
motions on a given scale is reduced by the same factor. 



4.6. Molecular Line Spectra toward GO.13-0.13 



Figure 8a-d show line emission from 16 molecular lines toward four positions in GO.13- 
0.13 corresponding to A to D in Table 1, respectively. The spectra of these pointed observa- 
tions show a rich chemistry similar to that seen throughout the CMZ. We focus only on five 
molecular lines HCN (1-0), HCO+(1-0), HNC(l-O), N 2 H+(l-0) and SiO (2-1), as Table 2 
show the peak velocities and peak intensities in T^ of all five spectral lines. Positions A and 
D fall in the E and NW boxes (see Fig. 5a). As a demonstration, the spectrum of these five 
molecular species toward position A is shown in Figure 9. The strong velocity component 
associated with GO.13-0.13 is centered around 60 km s" 1 (see Table 2). We also notice an 
additional weak velocity component between -50 km s _1 and km s _1 . 

Using RADEX0 (van der Tak et al. 20 7) 55 , we derived the column density of each 
molecular species to match the observed intensity with the assumption that molecular gas 
temperature and gas density are T= 200 K and 10 4 cm -3 , respectively and the velocity 
dispersion is 30 km s _1 . Table 3 shows the log of column density ratio of HCN/HNC, 
HCO + /N 2 H + and SiO/N 2 H + in columns 2-4 and the corresponding antenna temperature 
ratios in column 5-8 based on deep Mopra observations at the four positions A-D given 
column 1. Calibration uncertainties of 25% have not been included in tables 2 and 3. We 
note that position C lies at the center of the boot-shaped structure, cf. Fig. 2b, where the 
74 MHz emission peaks, implying that the cosmic ray ionization rate is higher than the other 
three positions. 



J http://www. sron.rug.nl/ vdtak/radex/index.shtml 



4.7. Chemical Modelling 



To examine the effect of high cosmic-ray ionization rates on chemistry, we make use of 
a time-dependent gas-grain chemical model, UCL_CHEM 31 to qualitatively investigate the 
behavior of the observed species. The UCL.CHEM is a gas-grain time-dependent model; 
for the purpose of this paper it is used as a two-stage model: Phase I follows the free 
fall collapse of a diffuse (10 2 cm -3 ) gas to a denser state (where the final density is a free 
parameter), while Phase II follows the chemistry as the gas and dust warm up due to either 
the increase of temperature and/or an enhanced cosmic ray ionization rate (Phase II). Apart 
from the very initial diffuse state (where the temperature is about 100K), during the collapse 
in Phase I the (coupled) gas and dust temperature remains constant at 10 K and atoms and 
molecules collide with, and freeze on to, grain surfaces. The advantage of this approach is 
that the ice composition is not assumed but it is derived by a time-dependent computation 
of the chemical evolution of the gas-dust interaction process which, in turns, depends on 
the density of the gas. Hydrogenation occurs rapidly on these surfaces, so that, for example, 
some percentage of carbon atoms accreting will rapidly become frozen out methane, CH 4 . 
The justification for a constant 10 K during most of Phase I is two-fold: first, in order for a 
diffuse gas to collapse under gravity and form a dense core, the temperature must remain low. 
Second, it is believed that methanol must form on icy mantles (and this implies an efficient 
freeze out of CO which can only occur for dust temperatures less than 20K) since experiments 
show that it can not form in the gas phase 56 (e.g., Geppert et al. 2006). In Phase II (which 
effectively represents the observed phase) the gas and dust temperature are decoupled; the 
temperature of the gas is varied from 50 to 200 K, while the dust temperature was kept to 20- 
30K. We employ the reaction rate data from the UMIST astrochemical database, augmenting 
it with grain-surface (hydrogenation) reactions 31,57 In both Phases nonthermal desorption is 
also considered 58 . 
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Figure 10 shows the abundance ratios of selected species as a function of time for a 
subset of our model grid where the final density was varied from 10 4 cm" 3 to 10 6 cm -3 
and the cosmic ray ionization rate from 1CT 15 s" 1 to 10~ 13 s _1 . The gas temperature in 
Phase II was assumed to be ~ 200K. Note, however, that while molecular ratios are in 
principle powerful tools in constraining the physical and chemical characteristics of a cloud, 
the assumption that all the selected species arise from the same region may well be incorrect. 
Nevertheless, we assume that the selected species are co-spatial and consider a comparison 
between the theoretical column density ratios with those derived from our observations (see 
Table 3). The [HCN/HNC] (black) and [HCO+/N 2 H+] (green) values, as drawn on the top 
right panel in dashed lines, agree with models of low density (10 4 cm~ 3 ) and high cosmic 
ray ionization rate ( = 1CT 13 s _1 H _1 . We note an increase in the gas density and or a 
decrease in the cosmic ray ionization rate does not change the abundance ratio [HCN/HNC] 
but increases the HCO + /N2H + (green) ratio to values well above the observed range of 0.005 
to 0.19 of Table 3. The cosmic ray ionization rate, inferred independently from this model, 
is remarkably similar to that estimated from Fel Ka line measurements (see Fig. 7a). In 
addition, the gas density predicted from this chemical model is consistent within the range 
of values inferred from LVG models, as shown in Figure 5b, c. 

It is known that Galactic center molecular clouds have high abundance of SiO 59 ' 60 ' 61 . 
In one estimate, Minh et al. (1992) find the SiO abundance relative to molecular hydrogen 
to be between ~ 10~ 7 and 10~ 8 . In other studies, Martin Pintado et al. (1997) also find 
a high abundance of SiO ~ 10~ 9 toward Galactic center molecular clouds. The observed 
[SiO/N2H + ] abundance ratio ranges between 0.27-0.44 and can not be matched by any 
model at late enough times that chemical equilibrium has been reached. The models predict 
even higher fractional abundance of [SiO/N 2 H + ]. In our models, we used an initial atomic 
abundance for Si depleted by a factor of 100 with respect to solar, as this represents the 
lower limit of measured initial abundances in a large sample of molecular clouds 62 . However, 
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all initial abudances 63 are then scaled to match a metallicity twice the solar value. Note that 
the abundance of SiO derived by our models is consistent with the fact that with high cosmic 
ray ionization rates the icy mantles are released and therefore Si is released in the gas phase. 
It is therefore puzzling to observe such a low fractional abundance in a gas where we expect 
cosmic ray ionization rates to be enhanced. There are three possibilities that could account 
for the low fractional abundance of [SiO/N2H + ] compared to that predicted by chemical 
model. One is that the medium is quite turbulent and hence chemical equilibrium is never 
reached and the gas is constantly being recycled; then the best solution must be found at 
earlier times (< 10 4 years). In this scenario all ratios can be qualitatively matched for the 
range of densities 10 4 " 5 cm -3 and the whole range of cosmic ray ionization rates investigated 
and no further constraints can be given by the models. The other is that the initial atomic 
abundance of Si is low. Although we already employ, in models shown in Figure 10, an 
abundance of the initial Si 1.64xl0~ 7 that is up to a factor of 100 depleted with respect to 
solar 64 . It could be that in these environments even more than a factor of 100 is depleted. 
Lastly, it is possible that our derivation of the column density from SiO (2-1) intensity is not 
accurate, e.g. optically thick. Thus, the intensity of high rotational transitions are needed 
to determine accurately the SiO column density. Additional detailed work is needed to sort 
out the discrepancy between the observed and modeled [SiO/N 2 H + ] abundance ratio. What 
is clear from our chemical modeling, as in shocks and and PDRs, a high SiO can also be 
generated by cosmic rays. 



4.8. Summary 

We showed that the molecular emission from a Galactic center cloud GO. 13-0. 13, a 
representattive molecular cloud in the CMZ, coincides with 74 MHz and Fel Ka line emission 
at 6.4 keV. This three-way correlation provides compelling evidence that relativistic electrons 
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are physically associated with the GO. 13-0. 13 molecular cloud. The high cosmic ray electron 
ionization rate ~ 10~ 13 s _1 H _1 is responsible for the Fel Ka line emission, heating the gas 
clouds to high temperatures and maintaining the velocity dispersion of clouds by increasing 
the wave damping time scale. One of the important conclusions of this study is that LVG 
modeling of multi-transition SiO observations gave molecular gas densities ~ 10 4-5 cm -3 and 
temperature ~ 100 — 1000K. The variation of the spectral index of nonthermal emission from 
GO. 13-0. 13 was shown to be unusually large Aa ~ 1.5. We suggested that the expansion 
of the GO. 13-0. 13 molecular cloud into nonthermal vertical filaments accelerates particles 
along the filaments 25 ' 48 , thus producing a young population of electrons with an unusually 
flat energy spectrum. 

Evidence for Fel Ka line and 74 MHz synchrotron emission from warm molecular gas 
supports the view that the Galactic center is a cosmic ray dominated region. We explored 
the chemistry in a region subject to high levels of cosmic ray ionization and showed that 
enhanced SiO emission, widespread throughout Galactic center molecular clouds, can be 
produced in the context of cosmic ray interaction with molecular gas. The model to explain 
enhanced SiO, as well as NH 3 and CH 3 OH, emission by cosmic ray driven chemistry is an 
alternative to the shock-driven chemistry throughout Galactic center molecular clouds. 
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Fig. 1.— (Top to Bottom) (a) A grayscale 74 MHz distribution of the inner 1.5° x 12' 
(Galactic coordinates I x b) convolved to a resolution of 2' x 2' based on VLA observations. 
The 74 MHz flux range is between -1 and 4 Jy. (b) Similar to (a) except the distribution of 
CO (3-2) integrated over ±200 km s _1 at 34" resolution 48 with flux range 100 and 2500 km 
s _1 K. (c) Contours of the CO (3-2) line emission shown in (b) are superimposed on the 74 
MHz image shown in (a), (d) Similar to (a) except observed with the GBT at 8.5 GHz with 
a resolution of 88" 34 . The flux range is between -5e-3 and 0.4 Jy. 



- 35 - 




Fig. 2. — (a - Left) Contours of CS (1-0) line emission integrated between and 50 km s _1 
from the GO. 13-0. 13 molecular cloud with a resolution of 45 //25 are superimposed on a 5 GHz 
continuum image with a spatial resolution of 10.8" x 5.5" (PA =-2. 6°). Contour levels are set 
at 2 to 10 km s" 1 K (T^) with 1 km s" 1 K interval, (b - Right) The distribution of CS (2-1) 
line emission (red) from GO. 13-0. 13 integrated over velocities between -12.1 < vlsr < 76.4 
km s -1 with a resolution of of 3.1" x 2.6" (PA=20.5°) is based on BIMA observations. The 
extended nonthermal filaments of the Arc at 1.4 GHz with a resolution of 10.7" x 10.1" 
(PA=28°) are shown in green. 
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Fig. 3.— (a - Top Left) Contours of 74 MHz emission with a resolution 122" x 64" (PA=-5°) 
with levels set at 0.8, 0.9, ...2.3 Jy beam -1 are superimposed on a grayscale image at 5 GHz. 
(b - Top Right) Contours of 327 MHz emission with levels (0.8, 0.9, ...2.3) x2 are shown 
with the same resolution as in (a) and are superimposed on the 5 GHz image. The white 
dashed lines show the location of cross cuts (see (d) below), (c - Left) The spectral index 
distribution a between 74 and 327 MHz emission. The 74 and 327 MHz are convolved to the 
same resolution and are background subtracted by 2 and 0.5 Jy, respectively. The bar at the 
bottom of the figure shows the color scale version of spectral index values ranging between 
-0.5 and 2 (see text), (d-e - Right) Top two panels show background subtracted intensity 
profiles at 327 and 74 MHz (top) and the corresponding spectral index distribution between 
these two frequencies (below) made from vertical latitude cross cut at constant longitude 
1=— 6', as drawn on (b) with a vertical dashed line. Bottom two panels are the same as the 
top two panels except that the cross cuts made horizontally at constant b=— 8'. Cross cuts 
are made on the 74 and 327 MHz images. 
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Fig. 3. — continued 
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Fig. 4. — Contours of Ka line emission at 6.4 keV line emission at (2, 3, 9)xl0 -6 ph s _1 
cm -2 beam -1 are superimposed on a CS (1-0) line emission, as shown in Figure 2b. CARMA 
image shows emission between 6.25 and 6.5 keV which is convolved to a Gaussian beam of 
30". 
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Fig. 5. — (a) Contours of SiO (5-4) line emission integrated between -20 and 100 km s~ 1 with 
levels at 8, 10, 12, 14, 16, 18, 20, 23 K km s _1 (T mb velocity) are superimposed on a 5 GHz 
grayscale image. Two rectangular boxes NW and E are drawn for LVG analysis, (b) The 
parameters of density njj cm~ 3 and temperature T (K) using LVG model are presented. 
SiO (5-4), (6-5), (2-1) line ratios are displayed in the bottom left corner in red, green and 
black, respectively. The dashed lines show the corresponding calibration uncertainties of the 
derived parameters at a level of 25%. The SiO line intensity ratios are extracted from the 
region presented by the E rectangular box, as drawn schematically on (a), (c) Similar to (b) 
except that data are from the NW box shown in (a), (d) Top three panels show the plots of 
line ratios as function of logarithm of pressure (nT) based on a grid of models. The observed 
line ratio values with their corresponding errors are drawn horizontally (green) in solid and 
dashed lines, respectively. The bottom panel shows the best x 2 fit to the grid of models. 
The vertical lines (green) are centered near minimum x 2 giving log(nT) values in the range 
of 6.4-6.9. 
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Fig. 5. — continued 
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Galactic Longitude 

Fig. 6. — Contours of SiO (2-1) line emission based on Mopra observations with levels 2, 
3.5, 5, 6.5, 8 K km s _1 are superimposed on a three color Spitzer/IRAC image combining 
Channels 4 (8/im) is red), 2 (4.5/im) is green and 1 (3.6/mi) is blue. The temperature scale 
is in T* A . The straight line in white shows the orientation of the nonthermal filament of the 
radio arc whereas the dashed circle marks roughly the size of the arc bubble. The positions 
where molecular abundances are measured are drawn as A, B, C and D. The box indicated 
the region mapped by Mopra. 
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Fig. 7. — (a - Top ) The variation of cosmic ray ionization rate as a function of magnetic field 
strength for different values of the spectral index of the radiation a. The black dot on each 
curve gives the value at which the magnetic field and particle energies are in equipartition. (b 
- Bottom) Solid curves show the total cooling rate for diffuse molecular gas for H 2 densities 
of 10 3 , 10 4 and 10 5 cm~ 3 . Rotational transitions of H2, CO and H2O, and by fine-structure 
transitions of 01 have been included; their individual contributions for 10 3 cm -3 are indicated 
by the dashed curves. The right hand axis shows the ionization rate by cosmic-ray electrons 
needed to supply the corresponding heating rate (see text). The vertical axis in both figures 
are logarithmic. 
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Fig. 8. — (a - Top Left 4) Deep pointed observations showing the spectra of 16 molecules 
toward position A in GO. 13-0. 13 (b - Top Right 4) Similar to (a) except for position B. (c - 
Top Left 4) Similar to (a) except for position C. (d - Top Right 4) Similar to (a) except for 
position D. 
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Table 1. Observed Positions within GO. 13-0. 13 



Source t b 
Name [°] [°] 



A 0.1007 -0.0792 

B 0.0632 -0.0709 

C 0.1132 -0.1084 

D 0.1299 -0.1292 
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Fig. 9. — Individual spectrum of five molecular species SiO, HCN, HNC, HCO + and N 2 H + 
toward position A based on Mopra observations. 
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Fig. 10.— Abundance ratios of HCN/HNC (black) and HCO+/N 2 H+ (green) and SiO/N 2 H+ 
(red) as a function of time in years for different values of cosmic ray ionization rate and 
molecular gas density. The chemical evolution of the GO. 13-0. 13 cloud uses the UCL.CHEM 
time-dependent gas-grain chemical model 31 . Using the same color scheme in the model, the 
observed abundance ratios for all four positions A-D are marked on the top right panel. 



-47- 



o tit 

i/3 ^ 



cc 

S-l 

O 

S-l 

H 



CD 
PI 

a 

cc 

cp 



o 
o 

I 

CO 

cc' 

o 



Pi 

o 



o 

Pj 

o 
> 

CP 
cc 

o 



CM 



O v; 



D 



O 

o 

I? 



O 

o 



O 



CO 

p 


CO 

p 


p 


CO 

p 


o 


d 


d 


d 


-H 


HH 


HH 


HH 


o> 


to 




X 


o 




CO 


CO 


















Ol 




cn 












HH 


HH 


HH 


o 

HH 






CO 




CN 


CD 


co 










CN 


ID 


5" 


CN 


■<# 


CN 






CO 




O 


O 




d 


6 


d 


d 


1 1 

TI 


1 1 
Ti 


1 1 


_u 

TI 


Ol 




Ol 


O 




CD 




o 


CN 


CN1 






o 


O 


o 


o 


Ol 


Ol 


Ol 


Ol 


o 


d 


d 




-H 


HH 


HH 


HH 


CO 




30 


o 




ID 


q 


00 


■X 1 


GO 




CO 


iO 




CO 


CO 










o 


o 


o 




o 


d 


d 




-H 


HH 


HH 


HH 


CN 


Ol 


CN 






ID 




Ol 


CN 


CN 




d 


O 


o 


O 


o 


Ol 


Ol 


Ol 


Ol 


o 


d 


d 


o 


-H 


-H 


HH 


-H 


CN 


CN 


ID 




Ol 




i—l 






d 


00 


d 


iO 


ID 


CN 


CO 




CO 


pn 




o 


o 


o 


Q 


o 


o 


o 


o 


+1 


HH 


HH 


HH 


Ol 


o 


CO 


X 




H 


CO 


0) 


CO 


CO 


CN 




o 


o 


o 


o 


Ol 




Ol 


Ol 


d 


d 


d 


d 


+1 


HH 


HH 


HH 




co 


CO 


a; 


Ol 


IO 


CO 


IO 




T— I 








IO 


CO 




CO 


CO 


CO 


CO 


p 


p 


p 


p 


d 


d 


d 


d 


HH 


HH 


HH 


HH 


cn 


a> 


IO 


X 


p 


°! 


CO 


OS 




CO 


CN 


I— I 



d 



o 

Cn 


o 

Ol 


o 

CC! 


d 


d 


d 


HH 


HH 


HH 


X 




IO 


X 


T— I 


CO 






CN 


IO 


CO 


CO 


CC 


O 


Q 



o 

CO 



o 


o 


CN 

o 


o 


o 


d 


d 


d 


+1 






HH 




to 




GO 




CO 


CN 


CO 


d 


d 


d 


d 



15! 



C~" 



o 

-H 



O 

o 



o 
o 

151 



■CO 

o 


CN 

o 


IO 

o 


o 


o 


d 


d 


d 


-M 


-H 






t- 


CO 


t- 




co 


lO 


CO 


q 


CN 




CN 


CN 



O 

in 



TO 


CN 

o 


CO 

o 


CO 

o 


o 


d 


d 


d 


+1 


4H 


HH 




o 


o 


H 






CN 


■* 




d 


d 


d 


d 



H3 

g 

i2f 



O 

o 
t-EL 

izf 



"H 
■CO' 



o 
o 

d 



35 



o 
d 

to 



O 
15 



o 
d 

-H 



" r-1 



o 
d 



q 
d 
HH 
us 

CN 

d 



3 

o 

CO 



<; m o q 



